FILLED JULIA SETS WITH EMPTY INTERIOR ARE COMPUTABLE 
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Abstract. We show that if a polynomial filled Julia set has empty interior, then it is 
computable. 



1. Introduction 

Julia sets of rational mappings. We recall the main definitions of complex dynamics 
relevant to our results only briefly; a good general reference is the book of Milnor |Milj . 
For a rational mapping R of degree degi? = d > 2 considered as a dynamical system on 
the Riemann sphere 

R : C -> C 

the Julia set is defined as the complement of the set where the dynamics is Lyapunov-stable: 

Definition 1.1. Denote Fr the set of points zeC having an open neighborhood U(z) on 
which the family of iterates R n \u^ z ) is equicontinuous. The set F R is called the Fatou set 
of R and its complement Jr = C \ Fr is the Julia set. 

In the case when the rational mapping is a polynomial 

P(z) = a + a x z + ■ • • + a d z d : C -> C 

an equivalent way of defining the Julia set is as follows. Obviously, there exists a neigh- 
borhood of oo on C on which the iterates of P uniformly converge to oo. Denoting A(oo) 
the maximal such domain of attraction of oo we have A(oo) C F(R). We then have 

J P = dA(oo). 

The bounded set C \ clv4(oo) is called the filled Julia set, and denoted K P \ it consists of 
points whose orbits under P remain bounded: 

K P = {ze t\ sup \P n (z)\ < oo}. 

n 

For future reference, let us list in a proposition below the main properties of Julia sets: 

Proposition 1.1. Let R : C — > C be a rational function. Then the following properties 
hold: 
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• Jr is a non-empty compact subset of<C which is completely invariant: R 1 (J(i?)) = 
J{R); 

• Jr = JRn for all n G N; 

• Jr is perfect; 

• if Jr has non-empty interior, then it is the whole of C; 

• let U C C be any open set with U fl Jr ^ 0. Then there exists n E N such that 
R n (U) D J R ; 

• periodic orbits of R are dense in Jr. 

Let us further comment on the last property. For a periodic point z = R p (zq) of period 
p its multiplier is the quantity A = X(z ) = DR p (z ). We may speak of the multiplier of a 
periodic cycle, as it is the same for all points in the cycle by the Chain Rule. In the case 
when | A] ^ 1, the dynamics in a sufficiently small neighborhood of the cycle is governed by 
the Mean Value Theorem: when |A| < 1, the cycle is attracting (super-attracting if A = 0), 
if |A| > 1 it is repelling. Both in the attracting and repelling cases, the dynamics can be 
locally linearized: 

(1.1) iP(R?(z))=A-iP(z) 

where ip is a conformal mapping of a small neighborhood of z to a disk around 0. By 
a classical result of Fatou, a rational mapping has at most finitely many non-repelling 
periodic orbits. The sharp bound on the number of such orbits is 2d — 2; it is established 
by Shishikura. Therefore, we may refine the last statement of Proposition 11.11 

• repelling periodic orbits are dense in Jr. 

In the case when |A| = 1, so that A = e 2md , 9 G K., the simplest to study is the parabolic 
case when 9 = n/m G Q, so A is a root of unity. In this case R p is not locally linearizable; 
it is not hard to see that Zq G Jr. In the complementary situation, two non-vacuous 
possibilities are considered: Cremer case, when RP is not linearizable, and Siegel case, 
when it is. In the latter case, the linearizing map ip from (jl.lj) conjugates the dynamics of 
RP on a neighborhood U(zq) to the irrational rotation by angle 9 (the rotation angle) on 
a disk around the origin. The maximal such rotation domain around zq is called a Siegel 
disk. 

A different kind of a rotation domain may occur for a non-polynomial rational mapping 
R. A Herman ring A is a conformal image 

v : {z G C| < r < \z\ < 1} ->• A, 

such that 

RPov(z) = v[e 2m6 z), 

for some p G N and 9 G R \ Q. 

To conclude our discussion of the basic facts of rational dynamics, let us formulate the 
Fatou-Sullivan Classification of the types of connected components of the Fatou set of a 
rational mapping. The term basin in what follows will describe the set of points whose 
orbits converge to a given attracting or parabolic periodic orbit under the iteration of R. 
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Fatou-Sullivan Classification of Fatou Components . For every connected component 
W C Fr there exists m G N such that the image H = R m (W) is periodic under the 
dynamics of R. Moreover, each periodic Fatou component H is of one of the following 
types: 

• a component of the basin of an attracting or a super- attracting periodic orbit; 

• a component of the basin of a parabolic periodic orbit; 

• a Siegel disk; 

• a Herman ring. 

Computability of subsets of IR n . The computability notions that we apply for subsets 
of M fc belong to the framework of computable analysis. Their roots can be traced to the 
pioneering work of Banach and Mazur of 1937 (sec [Maz ). The reader may find an extended 
exposition of the model of computation we are using in jBrCj . See [Wei] for a detailed 
discussion of the concepts of modern computable analysis. See also [BY] for a discussion 
of computability as applied to problems in Complex Dynamics. 

Given a compact set ScK*, our goal is to be able to approximate the set S with an 
arbitrarily high precision. Here we ought to specify what do "approximate" and "precision" 
mean in this setting. 

Denote B(x,r) C M fc the closed Euclidean ball with radius r centered at x. Denote C 
the set of finite unions US(^,r») C ~R k of closed balls whose radii and coordinates of the 
centers are all dyadic rationals: 



For a set K C M fc denote its e-neigborhood by B(K,e) = U S6 ^5(x, e). Recall that the 
Hausdorff metric on compact subsets of R fc is defined by 



That is, dn{X,Y) is the smallest quantity by which we have to "expand" X to cover Y 
and vice-versa. 

In computable analysis a compact set S C M. h is called computable, if there exists a 
Turing Machine (TM) M{m), such that on input m, M(m) outputs an encoding of C m G C 
satisfying dg(S, C m ) < 2~ m . That is, S is effectively approximable in the Hausdorff metric. 
This definition is quite robust. For example, it is equivalent to the computability of the 
distance function ds(x) = inf {\x— y\ : y G S}. In more practical terms, a set is computable 
if there exists a computer program which generates its image on a computer screen with 
one-pixel precision for an arbitrarily fine pixel size. 

In the case of computing Julia sets, we are not actually computing a single set, but a 
set- valued function that takes the coefficients of the rational function R(z) as an input, and 
outputs approximations of Jr. The machine computing Jr is given access to R(z) through 
an oracle that can provide an approximation of any coefficient of R(z) with an arbitrarily 




d H (X,Y) 



ini{e : X C B(Y,e) and Y C B(X,e)}. 
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high (but finite) requested precision. Such an oracle machine is denoted by . If P is 
the set of rational functions on which works correctly, then is said to be uniform 
on P. In case P is a singleton, the machine is said to be non-uniform. In other words, 
a uniform machine is designed to work on some family of inputs, whereas a non-uniform 
machine is only supposed to work on a specific input. 

We remark that Julia sets have appeared in a somewhat similar context in the book of 
Blum, Cucker, Shub, and Smale |BCSSj . They considered the question whether a Julia 
set is decidable in the the Blum-Shub-Smale (BSS) model of real computation. The BSS 
model is very different from the computable analysis model we use, and can be very roughly 
described as based on computation with infinite-precision real arithmetic. Some discussion 
of the differences between the models may be found in |BrC| and the references therein. It 
turns out (see Chapter 2.4 of |BCSSj ) that in the BSS model all, but the most trivial Julia 
sets are not decidable. More generally, sets with a fractional Hausdorff dimension, including 
the ones with very simple description, such as the Cantor set, are BSS-undecidable. 

Statements of Results. The main purpose of the present note is to prove the following 
result: 

Theorem 1.2. For each d > 2 there exists an oracle Turing machine , with the oracle 
representing the coefficients of a polynomial of degree d with the following property. Suppose 
that p(z) is a polynomial of degree d such that the filled Julia set K p has no interior (so 
that J p = K p ). Then K p is computable by M^. 

The known results in this direction, which make the statement of the theorem interest- 
ing, are the following. Independently, the second author |Brvj and Rettinger |Retj have 
demonstrated that hyperbolic Julia sets are polynomial-time computable by an oracle TM. 
In sharp contrast, in |BYj the second and third authors demonstrated the existence of 
non-computable Julia sets in the family of quadratic polynomials. The latter examples are 
given by polynomials with Siegel disks. It is shown in [BYJ that there exist non-computable 
examples with rather wild topology. By a method of Buff and Cheritat |BCj the boundary 
of the Siegel disk can be made smooth, and as a consequence, the critical point is not in 
the boundary and the Julia set is not locally connected (see |Roj for a discussion of the 
topological anomalies in such Julia sets). 

As our present theorem shows, however, the notions of topological complexity and com- 
putational complexity are rather distinct. It covers, for example, the case of Cremer 
quadratics, which are also known to possess non-locally connected Julia sets. Further, we 
find the following statement plausible: 

Conjecture. The procedure of [BY] can be carried out so as to yield non- computable and 
locally connected Siegel quadratics. 

Our main result implies that the reason one is not able to produce high resolution computer- 
generated images of Cremer Julia sets is due to the enormous computational time required 
(whether such images would be very informative is another question). In view of this, 
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it is of great interest to obtain a classification of Julia sets into classes of computational 
hardness. 

Finally, let us state a straightforward generalization of the main theorem to the case of a 
rational map without rotation domains (either Siegel disks or Herman rings). First note 
that by the Fatou-Sullivan Classification, every Fatou component of such a rational map 
belongs to the basin of a (super) attracting or parabolic orbit. A rational map / of degree d 
can have at most 2d — 2 of such orbits by Shishikura's result. Of these, the attracting orbits 
can be found algorithmically, with an arbitrary precision, as the condition \{f m )'{x)\ < 1 
can be verified with a certainty by a TM. For an algorithm to ascertain the position of a 
parabolic orbit a finite amount of information is needed. For example, we can ask for its 
period m, and a good rational approximation to distinguish the point among the roots of 
/"'(.:) z. 

Theorem 1.3. Let f be a rational map f : C — > C without rotation domains. Then its 
Julia set is computable in the spherical metric by an oracle Turing machine with the 
oracle representing the coefficients of f . The algorithm uses the nonuniform information 
on the number and positions of the parabolic orbits, as described above, as well as the 
(rational) values of the log's of the multipliers of the parabolic orbits. 

In particular, if there are no parabolic orbits the computation can be performed by a single 
oracle machine. As a function computed by an oracle Turing machine must be continuous 
on its domain (see e.g. JBrvj), we have the following corollary: 

Corollary 1.4. Denote 7Z d C C 2d ~ 2 the set of parameters of degree d rational maps without 
rotation domains or parabolic orbits. Then a Julia set depends continuously, in Hausdorff 
sense on C, on the parameter in lZ d . 

While this is well-known (see the discussion in |Doj ) . it is notable that the proof of this 
dynamical statement is produced by a computer science argument. 1 

Acknowledgement. The authors wish to thank John Milnor, whose encouragement and 
questions have inspired this work. 

2. Proof of the Main Theorem 

Idea of the proof. The idea of the argument is to find two sequences of Hausdorff bounds 
for the filled Julia set. The inside ones are given by rational approximation of the periodic 
cycles of p(z), accomplished by approximately solving a sequence of polynomial equations. 
The outside ones are given by preimages of a large disk centered at the origin. The algorithm 
stops when the difference between the two bounds is sufficiently small. This will always 
happen, provided the interior of the filled Julia set is empty. 

We make no claims about the efficiency of our algorithm. Indeed, in most cases it will 
be far too slow to be practical. The existence of an efficient algorithm for the challenging 
cases, such as Julia sets with Cremer points, is an interesting open problem. 

1 We owe this comment to J. Milnor 
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Lemma 2.1. For every natural n we can compute a sequence of rationals {qi} such that 

oo 

(2.1) B(J P , 2-^) m |J B( qi , 2'^) m B(J P , 2- n ). 

i=l 

Proof. First of all, note that for any polynomial Q(z), we can list rational approximations 
7*1, 7*2, . . . , r m of all the roots a>i, a 2 , ■ ■ ■ , a m of Q(z) — z with an arbitrarily good precision 
2~ k (a classical reference is Wey| ) . 

Let M > be some bound on |Q"(z)| in the area of the roots. Then |Q'( r ;)l > l+2~ k M 
implies that |Q'(a:j)| > 1, and in fact acts as a certificate that this is the case. It is easy to 
see that if, in fact, |Q'(ttj)| > 1, then eventually we will find a rational approximation 
of cti for which we know that |Q'(«i)| > 1. 

We now use the fact that the Julia set J p is the closure of the repelling periodic orbits 
of p. Thus we can compute a sequence {qi} of rational approximations of all the repelling 
periodic points such that: 

• For each i, there is a repelling periodic point cti with — a»| < 2~( n+3 \ and 

• for each repelling periodic point a, a 2~(™ +3 )-approximation of a will eventually 
appear in {qi}. 

We claim that (|2.1|) holds for the sequence {qi}. 

For any x G B(J P , 2~( n+2 )) there is a y G J p with \x — y\ < 2~( n+2 \ There is a repelling 
periodic point z with \z — y\ < 2~( n+3 \ and the rational point approximating it satisfies 
\z - qi\ < 2~( n+3 \ hence 

\x - fi | < \x - y\ + \y - z\ + \z - q { \ < 2" (n+2) + 2~ (n+3) + 2~ (n+3) = 2~ {n+1) . 

Soxe B{q u 2^ n+1 ^). 

On the other hand, for any there is a repelling z E J p with — z\ < 2~ (jl+3 \ hence 
B(q h 2"( n+1 )) C B(J p , 2"( n+1 ) + 2-("+ 3 )) C B(J P , 2~ n ). ^ □ 

Proof of Theorem ] 1. OA Fix m G N. Our algorithm to find C m G C works as follows. Evalu- 
ate a large number b > such that if \z\ > b then \p(z)\ > b and the orbit p n (z) —>■ oo. Set 
D = B(0, b). At the fc-th step: 

• compute the finite union B^ = U k =1 B(qi, 2 _ ( m+1 )) g C as in the previous lemma; 

• compute with precision 2~( m+3 ) the preimage p~ k (D), that is, find G C such 
that 

d H (TT k , P ~ Zk jD))<2-^ ] 

• if Dk C Bk output C m = Bk and terminate. Otherwise, go to step k + 1. 

We claim that this algorithm eventually terminates, and outputs an approximation C m G C 
with d H (C m ,K p ) < 2~ m . 

Termination. First, we note that p~ k (D) is compactly contained in p~^ k ~ l \D) for all 
k G N. Since in our case J p = K p , we have 

(2.2) K p — O p~ k {D) = {repelling periodic orbits of p}. 

fceN 
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By the previous lemma, there exists ko G N such that for all k > k$ we have 

B(K p ,2~^) m B k m B(K p ,2' m ). 

By (|2.2|) there exists k\ G N such that for all k > k\ we have 

p- k (D) C fi(K p ,2- (m+3) ), and hence D k C 2~ (m+2) ). 

Therefore, for > max(fco, &i) we have D k C and the algorithm will terminate. 

Correctness. Now suppose that the algorithm terminates on step k. Since D k C B k and 
K p C B(D k ,2~ (m+ V) we have K p C S(C m , 2~( m+3 )). On the other hand, U{ft} C K p , and 
thus 5 fe = C m C S(Kp, 2-( m+1 )). 

□ 

Proof of Theorem M.iA We will have to modify the way the set D k is generated in the 
preceding algorithm. Note that the procedure which finds periodic points of / of period 
m, and verifies the conditions \{f m )'{x)\ > 1, or |(/ m )'(a;)| < 1 with precision 2~ m will 
eventually certify every attracting periodic point. For each such point, we will also identify 
a Euclidean disk U x G C containing it with the property f m {U x ) <e U x . 

Further, for each parabolic point x of period I let Z7* G C be a topological disk compactly 
contained in the immediate basin of attraction of x (i.e. the union of attracting petals), 
with the properties: 

U* C U^ +l } and U fc U k x is an attracting petal. 

Such neighborhoods can easily be produced algorithmically given the non-uniform infor- 
mation on the parabolic orbits as stated in the Theorem. Our algorithm will then work 
as before, with the difference that the set D k is defined as a 2~( m+3 )-approximation of the 
complement of the preimage 

w k = r\ u u x )ur k ( u u k x ). 

x is attracting located so far x is parabolic 

By the Sullivan's Non- Wandering Theorem, we have J p = D(C \ W k ). The algorithm is 
verified exactly as in the previous theorem. □ 

3. Conclusion. 

We have shown that only rational maps with rotation domains (Siegel disks or Herman 
rings) may have uncomputable Julia sets. As was demonstrated earlier by the second and 
third authors |BYj . such a possibility is indeed realized. To understand the cause of this 
phenomenon, consider the practical problem of plotting the Julia set of a Siegel quadratic 
Po(z) = z 2 + e 2m6 z. The preimages of a large disk around the Julia set will rapidly converge 
to Kp g . However, when they are already quite close to Jp g , these preimages will still be 
quite far from the Siegel disk Ag, as they will not fill the narrow fjords leading towards 
dA e . To have a faithful picture of a 2~™-neighborhood of Jp g , one then needs to draw a 
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contour inside the Siegel disk, close to its boundary. In effect, this means knowing the 
inner radius of the Siegel disk. Indeed, we have the following easy statement: 

Proposition 3.1. Suppose the inner radius pg of A g is known. Then the Julia set J Pg is 
computable. 

Proof. The algorithm to produce the 2 _n approximation of the Julia set is the following. 
First, compute a large disk D around with Pe{D) s> D. Then, 

(I) compute a set Dk G C which is a 2~( n+3 )-approximation of the preimage P 9 ~ k (D); 
(II) set Wk to be the round disk with radius pg — 2~ k about the origin. Compute a set 

B k G C which is a 2~( n+3 )-approximation of P 9 ~ k (Wk); 
(III) if Dk is contained in a 2 - ( n+1 )-neighborhood of Bk, then output a 2~( n+1 )-neighborhood 
of D k — B k , and stop. If not, go to step (I). 

A proof of the validity of the algorithm follows along the same lines as that in Theorem ll.2l 

□ 

The well-recognizable images of the "good" Siegel Julia sets, such as the one for = 
(V5- l)/2, the golden mean, are indeed computable by this recipe: 

Proposition 3.2. The inner radius pg t is a computable real number. 

Proof. We appeal to the renormalization theory for golden- mean Siegel disks (see |McMj ) . 
which implies, in particular, that the boundary of Ag t is self-similar up to an exponentially 
small error. As a consequence, there exist C > 0, and A > 1 such that 

inf{|f(0)|, i = 0,...,q n }<pg t (l + C\- n ), 

where p n /q n is the n-th continued fraction convergent of 0*. □ 

In contrast, in the non-computable examples constructed in |BYj the boundary of the Siegel 
disk oscillates unpredictably at infinitely many scales. More precisely, recall the definition 
of the conformal radius r(U, u) of a simply-connected domain U with a marked point u 
in its interior. For the unique holomorphic bijection \l/ from the unit disk D to U with 
tf(0) = u, and *'(0) > set 

r(U,u) = $'(0). 

By the Koebe 1/4 Theorem of classical complex analysis, the inner radius of U as seen 
from u is at least \r{U, u). 

Now set rg = r(Ag,0). It is shown in |BYj that there exist parameter values r G R/Z 
for which the function 9 i— > rg is uncomputable on the set consisting of a single point {r}. 
On the other hand, we have: 

Proposition 3.3. Computability of Jp g implies computability of rg for Siegel parameter 
values 9. 

We note that a converse is also true, so computing the conformal radius is the only possible 
obstruction to drawing J Pg : 
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Figure 1. An image of the Siegel Julia set Jo* 
of Proposition 13.11 



produced by the algorithm 



Proposition 3.4. Suppose r$ is computable, then so is Jp g . 
We will require the following facts of complex analysis: 

Lemma 3.5. Let (U,u) be as above, and suppose U\ is a simply- connected subdomain of 
U containing the point u. Let e > be such that dU% C B(dU,e). Denote p and p\ the 
inner radii ofU and Ui respectively, as seen from u. There exist explicit positive constants 
Ci = Ci{r{U, u)), and C 2 = C 2 (r(U,u)) such that 

Ci\p-pi\ < \r(U,u)-r(Ux,u)\<C 2y fi. 

The second inequality is an elementary consequence of Koebe Theorem and the constant 
may be chosen as C 2 = A^r(U, u) (see |RZj . for example). The other inequality may be 
found in many advanced texts on Complex Analysis, e.g. [Pom]. 

Proof of Proposition We will show the computability of the inner radius of A$. The 
algorithm works as follows: 

(I) For k G N compute a set Dk G C which is a 2~ m -approximation of the preimage 

Pg k {D) ) for some sufficiently large disk D; 
(II) evaluate the conformal radius r(Dk,0) with precision 2~( m+1 ) (this can be done, 
for example, by using one of the numerous existing methods for computing the 
Riemann Mapping of a computable domain, see |Zhoj . for example), if necessairy, 
compute Dk or parts of D^ with a higher degree of precision; 
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(III) if r(Dk, 0) is 2 m GVclose to Tg then compute the inner radius p(Dk) around with 
precision 2~ m and output this number. Else, increment k and return to step (I). 

Termination. Let K = Kp be the filled Julia set of Pg. Then 

^ =0 D k = KDA e 
and Dq D Di D D% D . . .. Hence there will be a step k after which 

2~ 2m C? 

dist(dD k , J Pe ) < — — 2 — . 
dAg C Jp e , and by Lemma this implies that 

\r(D k , 0) - r(A e , 0)| = \r(D k , 0)-r g \< 2- m ^-C 2 = 2~ m C x , 
and the algorithm will terminate on step (III). 

Correctness. Now suppose the algorithm has terminated on step (III). Ag C D k , and 
Lemma fa. 51 implies that 

\p(D k )-p e \ <2~ m . 

□ 
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